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Abstract 

It is found that Lorenz systems can be unidirectionally coupled such that the chaos expands from 
the drive system. This is true if the response system is not chaotic, but admits a global attractor, 
an equilibrium or a cycle. The extension of sensitivity and period-doubling cascade are theoretically 
proved, and the appearance of cyclic chaos as well as intermittency in interconnected Lorenz systems are 
demonstrated. A possible connection of our results with the global weather unpredictability is provided. 
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1 Introduction 

In his famous study, to investigate the dynamics of the atmosphere, Lorenz |43| built a mathematical 
model consisting of a system of three differential equations in the following form, 

dxi 

—— = -axi + ax 2 
at 

dx2 -| \ 

=-XiXs + rxi - X 2 U-ij 

at 
dx3 

—— = XiX2 - 0X3, 

dt 

where cr, r and b are constants. 

System (ED is a simplihcation of a model, derived by Saltzman m, to study finite amplitude 
convection. The studies of Saltzman originate from the Rayleigh-Benard convection, which describes 
heat flow through a fluid, like air or water. In this modelling, one considers a fluid between two horizontal 
plates where the gravity is assumed to be in the downward direction and the temperature of the lower 
plate is maintained at a higher value than the temperature of the upper one. Rayleigh |55| found that if 
the temperature difference is kept at a constant value, then the system possesses a steady-state solution 
in which there is no motion and convection should take place if this solution becomes unstable. In 
other words, depending on the temperature difference between the plates, heat can be transferred by 
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conduction or by convection. Assuming variations in only xi — X 3 plane, Saltzman EH considered the 
equations 




d ['tp, 


d{xi,X3) 

d{ip,e) ATo dtp 


80 

dt ' d{xi,X 3 ) 


f)f) 

- ge- - vV^tp = 0 

UXl 


H dxi 


- = 0 , 


( 1 . 2 ) 


where '0 is a stream function for the two dimensional motion, 6 is the departure of temperature from that 
occurring in the state of no convection and the constants ATQ,H,g,e,v and n denote, respectively, the 
temperature contrast between the lower and upper boundaries of the fluid, the height of the fluid under 
consideration, the acceleration of gravity, the coefficient of thermal expansion, the kinematic viscosity 
and the thermal conductivity mini]- In his study, Saltzman EH achieved an infinite system by means 
of applying Fourier series methods to system (HH), and then used the simplification procedure proposed 
by Lorenz [32] to obtain a system with finite number of terms. Lorenz [33] set all but three Fourier 
coefficients equal to zero and as a consequence attained system uni), which describes an idealized model 
of a fluid. 

In system (11.11) . the variable X\ is proportional to the circulatory fluid flow velocity, while the variable 
xi is proportional to the temperature difference between the ascending and descending currents. Positive 
x\ values indicate clockwise rotations of the fluid and negative x\ values mean counterclockwise motions. 
The variable 0 : 3 , on the other hand, is proportional to the distortion of the vertical temperature profile 
from linearity, a positive value indicating that the strongest gradients occur near the boundaries. The 
parameters cr and r are called the Prandtl and Rayleigh numbers, respectively HHIllES]- 

The dynamics of the Lorenz system dm is very rich. For instance, with different values of the 
parameters cr, r and 5, the system can exhibit stable periodic orbits, homoclinic explosions, period¬ 
doubling bifurcations, and chaotic attractors |65| . 

The appearance of chaos in differential/discrete equations may be either endogenous or exogenous. 
As the first type of chaos birth, one can take into account the irregular motions that occur in Lorenz, 
Rossler, Chua systems, the logistic map. Duffing and Van der Pol equations [l8l IT^ [29l 1401 [58l 166] . To 
indicate the endogenous irregularity, we use: (i) ingredients of Devaney and Li-Yorke chaos, (ii) period¬ 
doubling route to chaos, (iii) intermittency, (iv) positive Lyapunov exponents. Symbolic dynamics and 
Smale horseshoes have been widely used for that purpose [121 EH 133 EZl SHI IMl EH ESj- While the 
endogenous chaos production is widespread and historically unique, the exogenous chaos as generated 
by irregular perturbations has not been intensively investigated yet. In our study, we will appeal to 
endogenous chaos, but mostly to exogenous chaos. 

In this paper, the main attention is given to the extension of chaos among interconnected Lorenz 
systems. We make use of unidirectionally coupled Lorenz systems such that the drive system is chaotic 
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and the response system possesses a stable equilibrium or a limit cycle. We theoretically prove that the 
chaos of the drive system makes the response system behave also chaotically. Extension of sensitivity and 
period-doubling cascade are rigorously approved. The appearance of cyclic irregular behavior is discussed, 
and it is shown that the phenomenon cannot be explained by means of generalized synchronization. 
Intermittency in coupled Lorenz systems is also demonstrated. 

The principal novelty of our investigation is that we create exogenous chaotic perturbations by means 
of the solutions of a chaotic Lorenz system, plug it into a regular Lorenz system, and find that chaos 
is inherited by the solutions of the latter. Such an approach has been widely used for differential 
equations before, but for regular disturbance functions. That is, it has been shown that an (almost) 
periodic perturbation function implies the existence of an (almost) periodic solution of the system. 
While the literature on chaos synchronization [Il[551[31[3Sl[3HlllZlElEn] has also produced methods 
of generating chaos in a system by plugging in terms that are chaotic, it relies on the asymptotic 
convergence between the chaotic exogenous terms and the solution of the response system for the proof 
of chaos creation. Instead, we provide a direct verification of the ingredients of chaos for the perturbed 
system m-m- Moreover, in Section M we represent the appearance of cyclic chaos, which cannot be 
reduced to generalized synchronization. Very interesting examples of applications of discrete dynamics to 
continuous chaos analysis were provided in the papers m-m- In these studies, the general technique 
of dynamical synthesis [H] was developed, and this technique was used in the paper g]. 

There are many published papers which have results about chaos considering first of all its mathe¬ 
matical meaning. This is true either for differential equations [njin] or data analysis [52]. Apparently 
there are still few articles with meteorological interpretation of chaos ingredients. We suppose that our 
rigorously approved idea for the extension of chaotic behavior from one Lorenz system to another will 
give a light for the justification of the erratic behavior observed in dynamical systems of meteorology. 

The question “Does the flap of a butterfly’s wings in Brazil set off a tornado in Texas ? " is very im¬ 
pressive and it has done a lot to popularize chaos for both mathematicians and non-mathematicians [45] . 
From this question one can immediately decide that the butterfly effect is a global phenomenon, and con¬ 
sequently, the underlying mathematics has to be investigated. Some of the authors say that the question 
relates sensitive dependence on initial conditions in dynamical systems considered as unpredictability 
for meteorological observations. Lorenz himself, in successive his talks and the book gS], was obsessed 
by the question and sincerely believed its possibility. He also supposed that his system can give a key 
for the positive answer of the question. Generally, analysis of chaotic dynamics in atmospheric models 
is rather numerical miMi [Ml sa uni or depend on the observation of time-series I2ai57|. In Section 
|8l we describe how one can use the rigorously approved results of the present paper to investigate the 
global behavior of the weather unpredictability. Our suggestions are not about a modelling, but rather 
an effort to answer the question why the weather is unpredictable at each point of the Earth, on the 
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basis of the Lorenz’s meteorological model and other models. We should recognize that all our discus¬ 
sions can be considered as a “toy object" in the theory, and according to the complexity phenomenon in 
meteorological investigations one can say that the investigation of chaos in meteorology still remains as 
a “toy object" [251115] . 


2 Coupled Lorenz Systems 


We couple Lorenz systems unidirectionally in such a way that the existing chaos propagates from one to 
another. We suppose that the coefficients cr, r and b are properly chosen in (HI so that the system is 
chaotic. In addition to system (HH), we consider another Lorenz system, 


dui 




dt 


= —aui + au 2 

= —U1U3 + rui — U2 
= U1U2 — bU 3 , 


(2.3) 


where the parameters a, r and b are such that the system is non-chaotic. That is, the system does not 
possess chaotic motions such that, for example, it admits a global asymptotically stable equilibrium or 
a globally attracting limit cycle. 

In order to realize the chaos extension, we perturb (12.311 with the solutions of (HH) to set up the 
system 

^ = -CT 2/1 +Wy 2 + gi{x{t)) 

-^ =-yiy3 + ryi-y2 +92ix{t)) (2.4) 

^ = yiJ/2 - bys + g3{x{t)), 

where x{t) = {xi(t),X 2 {t),X 3 {t)). The conditions on the continuous function g{x) = {gi{x), g 2 {x), g 3 {x)) 
is mentioned in the Appendix. 


3 Extension of Sensitivity 

In this section, we will demonstrate that the divergence of two initially nearby solutions (sensitivity) 
in the driving chaotic Lorenz system (HH) leads to the presence of the same feature in system (12a . 
Additionally, a third Lorenz system will be considered in order to show the maintainability of the process. 
The mathematical description of sensitivity and a theoretical proof for its extension are presented in the 
Appendix. 
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Let us take into account the system 


dt 

<m 

dt 

dys 


= —lOj/i + 10?/2 + 0.3xi{t) — 0.15 sin (a:i(t)) 
= -yiVs + 0.35yi - 2/2 + 1.6x2(t) 

= 2 / 12/2 - (8/3)2/3 + 0.1 tan(a;3(t)/65), 


(3.5) 


which is in the form of (12.411 with cr = 10, r = 0.35, b = 8/3, gi{x{t)) = 0.3xi(t) — 0.15sin(xi(t)), 
92 {x{t)) = 1. 6 x 2 ( 1 ) and 2 / 3 (x(l)) = 0.1 tan(x3(l)/65). Here, x{t) = (xi(l),X 2 (l),X 3 (l)) is a solution of 
(HI) with tr = 10, r = 28 and b = 8/3. The coefficients a, r and b are chosen such that chaos takes place 
in the dynamics of (HD gg. Besides, system (12.31) with the given values of cr, r, and b possesses a stable 
equilibrium point |65| . 

To reveal numerically the extension of sensitivity in system (13.51) . we represent in Figure [T] the 
projections of two initially nearby trajectories of the unidirectionally coupled system (HI T (13.51) on 
the yi — y2 — 2/3 space for t G [0,3]. The trajectory with blue color corresponds to the initial data 
xi(0) = —8.57, X2(0) = —2.39, X3(0) = 33.08, 2/1 (0) = 3.91, 2 / 2 ( 0 ) = 1.86, 2 / 3 ( 0 ) = 4.39, and the one with 
red color corresponds to the initial data xi(0) = —8.53, X2(0) = —2.47, X3(0) = 33.05, yi{0) = 3.91, 
2 / 2 ( 0 ) = 1.87, 2 / 3 ( 0 ) = 4.40. The divergence of the initially nearby trajectories seen in Figured) manifests 
the sensitivity feature in (13.51) . 



Figure 1: Extension of sensitivity in system (13.51) . The divergence of the initially nearby solutions of 
system (13.51) is observable in the figure. 


Now, we consider the system 


dt 

dz2 

dt 

dz3 


dt 


— lOzi + 102:2 + 122/1 (t) 

-^12:3 + O.lzi - Z2 + 20[2/2(t) + 2arctan(2/2(t)/5)] 
Z1Z2 - (8/3)z 3 - 82/3(t). 


(3.6) 


System (13.61) is also in the form of (12.41) . but this time the perturbations hi{y{t)) = \2yi{t)^ h 2 {y{t)) = 
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2 Q[y 2 {t) + 2 arctan(?/2(0/5)] and h 3 {y{t)) = — 8 y 3 {t) are provided by the solutions of (13.51) . 



Figure 2: Extension of sensitivity in system dSSl). 

Figure [2] shows the projections of two trajectories, which are initially nearby, of the 9—dimensional 
system (fTB+d^ + dXHl) on the zi — Z 2 — 23 space. The trajectory with blue color has the initial data 
a;i(0) = -8.57, 2 : 2 ( 0 ) = -2.39, 13 ( 0 ) = 33.08, ?/i(0) = 3.91, 2 / 2 ( 0 ) = 1.86, 2 / 3 ( 0 ) = 4.39, zi(0) = 6.92, 
^ 2 ( 0 ) = —6.18, 23 ( 0 ) = 10.48, whereas the one with red color has the initial data a;i(0) = —8.53, 
2 : 2 ( 0 ) = -2.47, 2 : 3 ( 0 ) = 33.05, 2 / 1 ( 0 ) = 3.91, 2 / 2 ( 0 ) = 1.87, 2 / 3 ( 0 ) = 4.40, 2 i( 0 ) = 6.89, 22 ( 0 ) = -6.18, 
23 ( 0 ) = 10.47. The utilized time interval is the same with Figure [TJ It is seen in Figure [5] that although 
the depicted trajectories are initially nearby, later they diverge from each other. In other words, it is 
demonstrated that the sensitivity of system (111311 is extended to dMD. Moreover, one can conclude from 
Figure [Hand Figured] that the system (I1.1I) + (I3.5|1 + (I3.6I) is also sensitive. 

In the next simulation, the trajectory of system (|ri|) + (I23])+(IS3|) with 2:1(0) = -12.89, 2:2(0) = 
-8.91, 2:3(0) = 36.59, 2/1(0) = -4.21, 2/2(0) = -4.96, 2/3(0) = 3.07, 2i(0) = -14.06, 22(0) = -8.38, 
■23(0) = 16.93 is considered. The three dimensional projections of the trajectory on the 2/1 ~ 2/2 ~ 2/3 
2i — 22 — 23 spaces are depicted in Figured] (a) and (b), respectively. Both of the pictures represented 
in Figure O manifest not only the chaos extension, but also the existence of a chaotic attractor in the 
9— dimensional phase space. It is worth noting that the projection on the 2:1 — 2:2 — 2:3 space is the 
classical Lorenz attractor [151155] . 

The next section is devoted to the extension of chaos obtained through period-doubling cascade. 


4 Extension of Period-Doubling Cascade 

Consider the Lorenz system (ED in which cr = 10, 6 = 8/3 and r is a parameter [51] |5S]. For the 
values of r between 99.98 and 100.795 the system possesses two symmetric stable periodic orbits such 
that one of them spirals round twice in 2:1 >0 and once in 2:1 < 0, whereas another spirals round twice 
in 2:1 < 0 and once in 2:1 > 0. The book |65j calls such periodic orbits as x^y and y^x, respectively. 
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a 


b 



Figure 3: The 3—dimensional projections of the chaotic trajectory of system (ll.lll + (l3.5|) + (l3.Hp . (a) The 
projection on the yi — y 2 — ys space, (b) The projection on the zi — Z 2 — Z 3 space. The picture represented 
in (a) illustrates the chaotic trajectory of the perturbed Lorenz system (I3.5L while the picture in (b) 
corresponds to the perturbed Lorenz system (13.61) . The pictures represented in (a) and (b) confirm both 
the extension of chaos and the existence of a chaotic attractor in the 9—dimensional phase space. 


More precisely, “x” is written every time when the orbit spirals round in Xi > 0, while “y” is written 
every time when it spirals round in xi < 0. As r decreases towards 99.98 a period-doubling bifurcation 
occurs in the system such that two new symmetric stable periodic orbits (x^yx'^y and y^xy^x) appear 
and the previous periodic orbits lose their stability El Eg. According to Franceschini m, system dnD 
undergoes infinitely many period-doubling bifurcations at the parameter values 99.547, 99.529, 99.5255 
and so on. The sequence of bifurcation parameter values accumulates at Too = 99.524. For values of r 
smaller than r^o infinitely many unstable periodic orbits take place in the dynamics of dTrli ElES]- 
To extend the period-doubling cascade of (H]), we take into account the system 


dyi 


di, 


= —lOyi -I- 10y2 + 1.8a;i(t) 


^ = -2/12/3 + 0.27yi - y2 + X2{t) 

at 

= 2/12/2 - (8/3)y3 -H 0.3x3(2), 


(4.7) 


where x{t) = (xi( 2 ), ^ 2 ( 2 ), ^ 3 ( 2 )) is a solution of ( 11 . 11 ) . Note that in the absence of driving, system (14.71) 
admits a stable equilibrium point, i.e., system (12.31) with a = 10, r = 0.27 and b = 8/3 does not admit 
chaos. 

By using Theorem 15.8 cni, one can verify that for each periodic x( 2 ), there exists a periodic solution 
of (14.71) with the same period. 

In Figure H we illustrate the stable periodic orbits of systems (11.11) and (14.7|) . Figure 01 (a) shows the 
y^x periodic orbit of (11.11) for r = 100.36, while Figured! (b) depicts the corresponding periodic orbit of 
system (14.71) . Similarly, Figure[ll (c) and (d) represent the y'^xy'^x periodic orbit of (11.11) with r = 99.74 
and the corresponding periodic orbit of (14.7L respectively. Figured] confirms that if (11.11) has a periodic 


7 




orbit, then (14.71) also has a periodic orbit with the same period. 






Figure 4: The stable periodic orbits of the systems dm and (liTD . 

Next, we continue with the extension of period-doubling cascade in Figure [5l The projection of the 
trajectory of system (11.11) with r = 99.51 corresponding to the initial data a;i(0) = 10.58, 0 : 2 ( 0 ) = 28.19, 
2 : 3 ( 0 ) = 53.32 on the 2:1 — X 3 plane is shown in Figure^ (a). Making use of the initial data j/i(0) = 2.23, 
2 / 2 ( 0 ) = 1.26, 2 / 3 ( 0 ) = 9.64, the projection of the corresponding trajectory of (14.71) on the 2/1 — 2/3 plane is 
depicted in Figure [5l (b). Moreover, the irregular behavior of the 2/3 coordinate over time is illustrated 
in Figure [HI The simulation results reveal that the period-doubling cascade of (11.11) is extended to (14.71) . 
A theoretical investigation of the extension of period-doubling cascade is provided in the Appendix. 

a b 




Figure 5: Extension of period-doubling cascade in the unidirectionally coupled Lorenz systems ini) + 

dm- 


In the next section, the extension of intermittency is considered. 
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Figure 6 : The irregular behavior of the 2/3 coordinate of system (I4.7I) with r = 99.51. 


5 Extension of Intermittency 


Pomeau and Manneville [S3] observed intermittency in the Lorenz system ( 11 . 11 ) . where < 7 = 10 ,5 = 8/3 
and r is slightly larger than the critical value Tc ~ 166.06. Let us use r = 166.3 in system (ED such that 
intermittency is present. We perturb system (12.3L where a = 10, r = 0.35, b = 8/3, with solutions of 
(HD, and set up the following system. 


^ = -lOj/i 
dt ^ 

dh 


10 y2 + 0.7xi{t) 


= -2/12/3 + 0.35?/i -y2- 2^2 (0 


dh 


2/3 

dt 


= 2/12/2 - ( 8 / 3 ) 2/3 + 0.2x3(2). 


(5.8) 


The graphs of the yi, 1/2 and 1/3 coordinates of (15.81) are shown in Figure[3by making use of the initial 
data xi(0) = —23.3, X2(0) = 38.3, ^ 3 ( 0 ) = 193.4, 2/i(0) = 1.3, 2 / 2 ( 0 ) = 5.5, 2 / 3 ( 0 ) = 12.1. It is revealed in 
Figure [7| that regular oscillations are interrupted by irregular ones, i.e., the intermittent behavior of the 
prior Lorenz system is extended even if system (12.31) admits a stable equilibrium point. 



t 


Figure 7: Extension of intermittency in system (15.81) . The behaviors of the 2 / 1 , 2/2 and 2/3 coordinates are 
shown in pictures (a), (b) and (c), respectively. The extension of the intermittent behavior is observable 
such that regular motions are interrupted by irregular ones. 

For some parameter values, the Lorenz system can exhibit limit cycles m- In the next section, we 
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will consider the Lorenz system (12.31) with a globally attracting limit cycle and verify numerically how 
to achieve a motion that behaves chaotically and cyclically in the same time. 


6 Cyclic Chaos 


In our previous illustrations, we considered system (lOll with a stable equilibrium point. Now, we consider 
the model with a limit cycle. The numerical simulations represented in this section are theoretically 
based on the paper |10j . where the main result is about the existence of infinitely many unstable periodic 
solutions and extension of sensitivity. 

Let us consider the systems (ED and with the coefficients cr — 10, v — 28, h — 8/3 and a — 10, 
r = 350, b = 8/3, respectively, such that (11.111 is chaotic and (12.31) possesses a globally attracting limit 
cycle m- We perturb system (12.311 with the solutions of (II.IL and constitute the system 


^ = -lOm 

dt ^ 

ay 


102/2 + 2.3a:i(<) 




ay 


dt 


= -yiVs + 3502/1 - 2/2 + 2x2(1) 
= 2/12/2 - (8/3)2/3 + 1.5x3(1). 


(6.9) 


Making use of the solution of (ED corresponding to the initial data xi(0) = 5.71, X2(0) = 9.01, 
X3(0) = 17.06, we depict the trajectory of (16.911 with 2/i(0) = —21.67, 2 / 2 ( 0 ) = 34.33, 2 / 3 ( 0 ) = 346.38 in 
Figure [HI (a). The projection of the same trajectory on the 2/1 — 2/2 plane is shown in Figure |S1 (b). It is 
seen in both figures that the trajectory behaves chaotically around the limit cycle of (ED- 

a b 


150 



Figure 8: The projections of the chaotic trajectory produced by the coupled system (II.ID + (I6.9I) . (a) The 
3—dimensional projection on the 2/1 ~ 2/2 ~ 2/3 space; (b) The 2—dimensional projection on the 2/1 — 2/2 
plane. The pictures in (a) and (b) represent a motion that behaves both chaotically and cyclically around 
the stable limit cycle of (12.31) . 


To confirm one more time that the trajectory considered in Figure [5] is essentially chaotic, the graph 
of the 2/2 coordinate of system (16.91) is illustrated in Figure [HI Although system (12.31) possesses a globally 
attracting limit cycle, the simulations seen in Figure|5]and FigurelHlindicate that the applied perturbation 
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makes the system (16.91) behave chaotically. In other words, the chaotic behavior is seized by the limit 
cycle of system (lOl) . and as a result a motion which behaves both chaotically and cyclically appears. 
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Figure 9: The time-series for the ?/ 2 —coordinate of system f|6.9[) . The picture confirms the chaotic behavior 
of the motion. The remaining coordinates of system (|6.9p . which are not just pictured here, behave also 
chaotically. 


In order to compare our approach with that of generalized synchronization [UllSllMlIMlEni, let US 
apply the auxiliary system approach [1] [53] to the couple (ITT]) + (lOl) . 

The corresponding auxiliary system is 


dz 

= —lOzi + 1002 + 2.3xi(t) 
at 
dz2 


=-Z 1 Z 3 + 350zi - + 2a;2(t) 

^ = ^ 1-22 - ( 8 / 3 ) 2;3 + 1.5X3{t). 
at 


i 


( 6 . 10 ) 


The projection of the stroboscopic plot of the 9—dimensional system (11.11) + (16.91) + (16.101) on the 
2/2 — Z 2 plane is depicted in Figure 1101 The figure is obtained by marking the trajectory with the 
initial data a:i(0) = 5.71, a; 2 ( 0 ) = 9.01, ^ 3 ( 0 ) = 17.06, 2/1 (0) = —21.67, 2 / 2 ( 0 ) = 34.33, 2 / 3 ( 0 ) = 346.38, 
zi(0) = —46.26, 2 : 2 ( 0 ) = —49.73, 23 ( 0 ) = 415.87, and by omitting the first 200 iterations. It is observable 
in Figure[Tn]that the stroboscopic plot is not on the line Z 2 = 2 / 2 - Therefore, we conclude that generalized 
synchronization does not take place in the dynamics of the couple (fTTl) + (lOl). 



^2 


Figure 10: Application of the auxiliary system approach to system (11.11) + (16.91) indicates that generalized 
synchronization does not exist for the couple. 


Another approach to investigate the presence or absence of generalized synchronization is the evalu- 
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ation of conditional Lyapunov exponents |25l 1551155] . 

To determine the conditional Lyapunov exponents, we take into account the following variational 
equations for system (16.91) . 




dt 

dt 

*3 


-106 +106 

(-2/3(6 + 350)6 - 6 - 2/1(66 
2/2(66 + 2/1(66 - ( 8 / 3 ) 6 - 


( 6 . 11 ) 


Utilizing the solution y{t) of (16.91) corresponding to the initial data xi(0) = 5.71, ^ 2 ( 0 ) = 9.01, X 3 ( 0 ) = 
17.06, yi(0) = —21.67, 1 / 2 ( 0 ) = 34.33, 1 / 3 ( 0 ) = 346.38, the largest Lyapunov exponent of system (16.111) is 
evaluated as 0.0226. That is, system (16.91) possesses a positive conditional Lyapunov exponent, and this 
result reveals one more time the absence of generalized synchronization. 


7 Self-organization and Synergetics in Lorenz Systems 


To illustrate the extension of chaos in large collections of interconnected Lorenz systems, let us introduce 
the following 27—dimensional system consisting of the subsystems Si^ S 2 , ■ ■ ■, Sq : 


dxi 

dt 

dX 2 

dt 

dX 3 

dt 


— lOxi + 10 a ;2 


—X 1 X 3 + 28a;i — X 2 




3^12:2 - (8/3)a;3 


dyi _ 
dt 

dy 2 _ 

dt 

dys _ 
dt 

dz\ 

dt 

dZ 2 

dt 

dZ 3 

dt 

dwi 

dt 

dW 2 

dt 

dW 3 

dt 


-lOyi + IO 1/2 + 8 xi{t) 

-2/12/3 + O.2I1/1 - 1/2 + X 2 it) + 0 . 001 x 2(6 
1/12/2 - (8/3)i/3 + 2 x 3(6 


>^2 


= -lOzi + 10Z2 + 4X2 (t) 

= —Z 1 Z 3 + 0 . 02 zi — Z2 + 3x3 (t) 

= Z1Z2 - (8/3)z 3 + xi(6 + 0.1 cos(xi(6) 




= — lOlCi + 10 w2 — 

= —W1W3 + 0 . 34 ii;i 


Xi(t) 

— 1 C 2 + 4 tanh(x 2 ( 6 ) 


5*4 


= W1W2 - { 8 / 3 )w 3 - 5x3(6 
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dt 

dt 

d-Cs 

dt 

drji 

dt 

dm 

dt 

dm 

dt 

dKi 

dt 

dK2 

dt 

dm 

dt 

dpi 

dt 

dp2 

dt 

dp3 

dt 


= -lOCi + IOC 2 + tan(yi(t)/20) 
= -CiC3 + 0.12Ci-C2 +2.52/2(0 
= C1C2 - ( 8 / 3 )C 3 - 102/3(0 

= — lOr/i + 107/2 + 82/1(0 

= -mm + 0.297/1 - 7/2 + 4 . 52 / 3(0 

= Vim - (8/3)7/3 - e*'^0)/30 

= —IOki + 10k 2 + 4zi(0 
= —Ki7t3 + 0.19 ki — K2 + 9^2(0 
= mm - (8/3)k3 + 623(0 

= -10/71 + 10p2 + 47771 (0 
= -PlP 3 + 0.17pi - P2 + 77772(0 
= Pi/72 - (8/3)p 3 - 3tanh(7773(i)) 






S7 


^^8 


= — lOri + 10r2 + arctan( 777 i(t)) 
^ = -rir3 + 0.32ri - T2 + 9'W2it) 
^ = T1T2 - ( 8 / 3 )r 3 + 7773(0. 


>^9 


The coefficients of Si are chosen in such a way that the system is chaotic [33]. The systems 
S2, S3,..., Sg are designed such that if the corresponding perturbations x{t), y(t), z{t), w(t) are chaotic, 
then the systems possess chaos. However, in the absence of the perturbations, S2, S3,..., Sg admit 
stable equilibria and they are all non-chaotic. The connection topology of the systems Si, S2, ■ ■ ■, Sg is 
represented in Figure [11] On the other hand, Figure |T3] depicts the chaotic attractors corresponding the 
each Si, i = 1,2,...,9, such that collectively the picture can be considered as the chaotic attractor of 
the whole 27—dimensional system. One can confirm that Figure [12] supports our ideas such that the 
chaos of Si generates chaos in the remaining subsystems even if they are non-chaotic in the absence of 
the perturbations. 

The idea of the transition of chaos from one system to another as well as the arrangement of chaos in an 
ordered way can be considered as another level of self-organization [sniEa. Durrenmatt m described 
that “... a system is self-organizing if it acquires a spatial, temporal or functional structure without 
specific interference from the outside. By “specific" we mean that the structure of functioning is not 
impressed on the system, but the system is acted upon from the outside in a nonspecific fashion." There 
are three approaches to self-organization, namely thermodynamic (dissipative structures), synergetic and 
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Figure 11: The topology of the chaos extension through the interconnected Lorenz systems Si, S2, ■ ■ ■, Sg. 


the autowave. For the theory of dynamical systems (e.g. differential equations) the phenomenon means 
that an autonomous system of equations admits a regular and stable motion (periodic, quasiperiodic, 
almost periodic). In the literature, this is called as autowaves processes [T^] or self-excited oscillations 
m- We are inclined to add to the list another phenomenon, which is a consequence of the chaos 
extension. Consider the subsystems Si, S2, ■ ■ ■, Sg once again (in general, arbitrary finite number of 
systems can be considered). Because of the connections and the conditions discovered in our analysis 
(see Appendix), likewise Si all the other subsystems. Si, i = 2,3,... ,9, are also chaotic. We suppose that 
this is a self-organization. This phenomenon can be restricted only for autonomous systems or it can be 
even interpreted for non-autonomous systems, too. So, we can say that the extension of unpredictability 
is an example of self-organization, that is a coherent behavior of a large number of systems m- 

In his fascinating paper, the German theoretical physicist Hermann Haken m introduced a new inter¬ 
disciplinary field of science, synergetics, which deals with the origins and the evolution of spatio-temporal 
structures. The profound part of synergetics is based on the dynamical systems theory. Depending on 
the discussion of our manuscript, it is natural that we concentrate on the differential equations, and 
everything that will be mentioned below about synergetics concerns first of all dynamical systems with 
mathematical approach. One of the main features of systems in synergetics is self-organization, which 
has been discussed above, and we approved that the phenomenon is present in the extension of chaos 
among Lorenz systems. According to Haken m, the central question in synergetics is whether there 
are general principles that govern the self-organized formation of structures and/or functions. The main 
principles by the founder of the theory are instability, order parameters, and slaving |30| . Instability is 
understood as the formation or collapse of structures (patterns). This is very common in fluid dynamics, 
lasers, chemistry and biology 1301 [311 ISlISIlill- A number of examples of instability can be found in the 
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literature about morphogenesis |67j and the pattern formation examples can be found in fluid dynamics. 
The phenomenon is called as instability because the former state of fluid transforms to a new one, loses 
its ability to persist, and becomes unstable. We see instability in the chaos extension, as consecutive 
chaotification of systems S 2 , S 3 ,..., Sg joined to the source Si of chaos. The concepts of the order 
parameter and slaving are strongly connected in synergetics. For differential equations theory, order 
parameters mean those phase variables, whose behavior formate the main properties of a macroscopic 
structure, which dominate over all other variables in the formation such that they can even depend on 
the order parameters functionally. The dependence that is proved (discovered) mathematically is what 
we call as slaving. It is not difficult to see in the chaos extension mechanism that the variables of the 
system are order parameters, and they determine the chaotic behavior of the joined systems’ variables. 
That is, the slaving principle is present here. 

The next section is devoted to the possible connections of our results about interconnected Lorenz 
systems with the global weather unpredictability. 


8 Connection with Global Weather Unpredictability 

Lorenz |43| was the first who discovered sensitivity with the aid of system (El) and then made the 
conclusion on the butterfly effect. Nowadays, there is an agreement that the butterfly effect exists, if 
we mean sensitivity=unpredictability. It seems that Lorenz himself believed that sensitivity discovered 
in his equation is a strong indicator of the butterfly effect in its original meteorological sense. Possibly 
his intuition is based on the idea that the system of ordinary differential equations is derived from a 
system of partial differential equations. There should be a deeper interpretation for the effect of chaotic 
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dynamics in the three dimensional system on the infinite dimensional one. We also believe that the 
opinion of Lorenz, who considered his results as an evidence of the meteorological butterfly effect, is very 
reasonable. Moreover, his claim has to be considered as a challenging problem for mathematicians. If 
one thinks positively on the subject, then by our opinion several next questions emerge. The first one is 
whether sensitivity in meteorological models is a reflection of the butterfly effect. Definitely, this question 
needs a thoroughly investigation. Possibly it requests a deep analysis on the basis of ordinary and partial 
differential equations. The problem is not solved in this section at all. We axiomatize somehow the state 
assuming that the butterfly effect is sensitivity in the mathematical sense or, more generally, chaos. 
We can also reduce the question by considering the problem of unpredictability through sensitivity. 
Consequently, the following questions are reasonable: Can one explain the global unpredictability of 
weather by applying models similar to the Lorenz system? How can Lorenz systems be utilized for a 
global description of the weather? 

The physical properties of the atmosphere are not the same throughout the Earth. The tropical 
atmosphere possesses considerably different behavior from those in the temperate and polar latitudes, 
as if it were a different fluid m- Taking inspiration from its multifaceted structure, we propose to 
consider the atmosphere divided into subregions such that the dynamic properties of each region differs 
considerably from the others. In this case, one can suppose that the dynamics of each subregion of the 
atmosphere subjects to its own Lorenz system. That is, for different subregions the coefficients of the 
corresponding Lorenz system are different. Since for some parameter values chaos can take place in the 
Lorenz system and for some not, such chaotic or non-chaotic motions should have prolonged forever, 
conflicting the realistic dynamics of the atmosphere, where global unpredictability is present. To extend 
our attitude for the butterfly effect, we propose that instability, which may occur in a subregion, can be 
imported to neighbor subregions of the atmosphere, such that chaos occurs not only endogenously, but 
also exogenously. In other words, exterior perturbations influencing a part of the atmosphere may cause 
a chaotic behavior to occur in that region. In addition to this, we suppose that these perturbations most 
probably originate through the neighboring regions within the atmosphere, and the dynamics of coupled 
Lorenz systems can help to analyze this. 

The results presented in the previous sections can be useful to investigate the underlying reasons of 
the global weather unpredictability under the following assumptions: 

(i) The whole atmosphere of the Earth is partitioned in a finite number of subregions; 

(ii) In each of the subregions the dynamics of the weather is governed by the Lorenz system with 
certain coefficients; 

(iii) There are subregions for which the corresponding Lorenz systems admit a chaos with the main 
ingredient as sensitivity, which means unpredictability of weather in the meteorological sense, and 
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there are subregions, where Lorenz systems are non-chaotic and with equilibriums or cycles as 
global attractors; 

(iv) The Lorenz systems are connected unidirectionally. 

Let us localize the global process by taking into account only two adjacent subregions of the atmo¬ 
sphere, labeled A and B. In the beginning, the subregion A is assumed to be chaotic, while the subregion 
B is non-chaotic. By the phrase “chaotic subregion" we mean that the coefficients of the corresponding 
Lorenz system are such that the system possesses a chaotic attractor. In a similar way, one should un¬ 
derstand from the phrase “non-chaotic subregion" that the corresponding Lorenz system does not exhibit 
chaotic motions such that, for example, it admits a global asymptotically stable equilibrium or a globally 
attracting limit cycle. 

Suppose that the dynamics of A is described by the chaotic Lorenz system (inj. Besides, system (El 
represents the dynamics of B after the transmission of chaos, and system (El represents the dynamics 
of B before the process is carried out. According to the theoretical results of the present paper, the 
chaos of the subregion A influences the subregion B in such a way that the latter also becomes chaotic 
even if it is initially non-chaotic. The propagation mechanism is represented schematically in Figure [H 
By chaos propagation, we mean the process of unidirectional coupling of Lorenz systems. Figure [H (a) 
illustrates the dynamics during the transmission of chaos. After the transmission of unpredictability is 
achieved, the dynamics of both subregions, A and B, exhibit chaotic behavior as shown in Figure [H 
(b). 


Chaotic 

Propagation 
of Chaos 

Non-chaotic 

Subregion A 


Subregion B 


b 



Figure 13: Schematic representation of the chaos extension mechanism, (a) The dynamics during the 
transmission of chaos, (b) The state of the weather after the transmission of chaos. The unidirectional 
coupling of Lorenz systems gives rise to the chaotification of the initially non-chaotic system such that 
as a result the unpredictability has been propagated from subregion A to subregion B. 

The mentioned local process can be maintained by considering more subregions, whose dynamics are 
also described by Lorenz systems. For instance, one can suppose that the systems Si, 82 , ■ ■ ■, Sg defined 
in Section [7] reflect the dynamics of the weather in nine different subregions. The applied perturbations 
may not be in accordance with realistic air flows in the atmosphere. However, the exemplification 
reveals the propagation of unpredictability and indicate the possibility for the usage of different types of 
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perturbations in the systems. We do not take into account changes which may happen because of the 
day light evolution, variety of land forms, seasonal differences in the region etc., but what we propose 
is to connect regional mathematical models into a global net so that understanding the unpredictability 
becomes possible. We make use of “toy" perturbations due to the lack of preexisting ones, which should 
be found through experimental investigations. What we propose in this section is a small step in the 
mathematical approach to the complexity of the weather. 

It is worth noting that the chaotification principles proposed in this paper are not specific for the 
Lorenz system (see Appendix), and they can be applied to other meteorological models as well, without 
any restrictions on the dimension and the number of the coupled systems. For example, one can consider 
the Lorenz model of general circulation of the atmosphere |44) 


^ = -X^ - - dXi + dF 

^=XiX2-bXiX3-X2 + G ( 8 . 12 ) 

at 

= bXiX2 + X1X3 — A3, 
at 


where Xi represents the strength of a large scale westerly wind current, X 2 and X 3 represent the cosine 
and sine phases of a chain of superposed large-scale eddies, the parameter F represents the external¬ 
heating contrast, and G represents the heating contrast between oceans and continents. The coefficient 
b, if greater than unity, allows the displacement to occur more rapidly than the amplification, and 
the coefficient a, if less than unity, allows the westerly current to damp less rapidly than the eddies 


For a > 0 and b > —1, let us take into account the Lyapunov function V{X) = X^ -|- A| -|- A|, 
and set a = min {a, 1} , 6 = -t- G^. One can verify that < ~2a (Aj -|- A| -|- A|) -|- 

2by/Xf + A| + A| and = 2^JXf + A| + A|. Therefore, the conditions of Theorem llO.il which 


dX 


is mentioned in the Appendix, are satisfied with a(r) = r^, b{r) = 2ar^ — 2br, c{r) = 2r, Mq = 1 and 

1 + 6 


B = 


Consequently, our theory is also applicable to the Lorenz model of general circulation of the 


atmosphere. 


9 Conclusion 

In the present study, we investigate the dynamics of unidirectionally coupled Lorenz systems. It is 
rigorously proved that chaos can be extended from one Lorenz system to another. The extension of 
period-doubling cascade and sensitivity, which is the main ingredient of chaos, are shown both theo¬ 
retically and numerically. Besides, the emergence of cyclic chaos, intermittency, and the concepts of 
self-organization and synergetics are considered for interconnected Lorenz systems. The results are valid 
if the drive Lorenz system is chaotic and the response system is non-chaotic, but admits a global asymp- 
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totically stable equilibrium or a globally attracting limit cycle. Our approach can give a light on the 
question about how global weather processes have to be described through mathematical models. A 
possible connection of the presented results with the global weather unpredictability is provided in the 
paper. The usage of our approach for the investigation of global weather unpredictability is a one more 
small step in the mathematical approach to the complexity of the weather. This is not a modelling of the 
atmosphere, but rather an effort to explain how the weather unpredictability can be arranged over the 
Earth on the basis of the Lorenz’s meteorological model. In fact, this is also true for other meteorological 
models, since mathematical properties of stability, attraction and chaotic attractors are common for all 
models. It is shown that our results can be used for the Lorenz model of general circulation of the atmo¬ 
sphere |44| too. The question whether the overlapping of two chaotic dynamics may produce regularity 
can also be considered in future investigations. We guess that it is not possible, but an analysis has to 
be made. 
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10 Appendix: The Mathematical Background 

In our theoretical discussions, we consider more general coupled systems, which are not necessarily Lorenz 
systems. We will denote by R and N the sets of real numbers and natural numbers, respectively, and we 
will make use of the usual Euclidean norm for vectors. 

Let us consider the autonomous systems 


dx 

dt 


F{x), 


(10.13) 


and 


du 

dt 


fiu), 


(10.14) 


where t > 0 and the functions F : —>■ R™ and / : R" —>■ R” are continuous in their arguments. 
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We perturb system (110.141) with the solutions of (110.131) and obtain the system in the form, 


-^ = f{y) + fj-gixit)), (10.15) 

where the real number /i is nonzero and the function g : R"* —> R" is continuous. It is worth noting that 
the systems (II.IL (12.31) and (12.41) are in the form of (|10.13l) . (110.141) and (110.151) . respectively. 

We mainly assume that system (110.131) possesses a chaotic attractor, let us say a set in R"*. Fix xq 
from the attractor and take a solution x{t) of (110.131) with x(0) = xq. Since we use the solution x{t) as 
a perturbation in (I10.15|) . we call it as chaotic function. Chaotic functions may be irregular as well as 
regular (periodic and unstable) [3TJ 031021E31 EH ES] • 

Our purpose is the prove rigorously the extension of chaos from system (110.131) to system (110.151) . In 
our theoretical discussions, we request the existence of a bounded positively invariant region for system 
(110.151) . Such an invariant region can be achieved by different methods and one of them is mentioned 
in the next part. We will show the extension of sensitivity and the existence infinitely many unstable 
periodic solutions in Subsections 110.21 and 110.31 respectively. 

In the following parts, for a given solution x{t) of system (I10.13p . we will denote by ^O) Vo) the 

unique solution of system (110.151) satisfying the initial condition <(>a:(t) (toj to, yo) = Ho- 

10.1 Existence of a bounded positively invariant region 

Making benefit of Lyapunov functions and uniform ultimate boundedness [MUZO], we present a method 
in Theorem 1 10.1 1 for the existence of a bounded positively invariant set of system (110.151) . Then, we will 
apply this technique to the Lorenz system. 

Solutions of system (I10.15|) are uniformly ultimately bounded if there exists a number Bq > 0 and 
corresponding to any number a > 0 there exists a number T{a) > 0 such that ||yo|| < a implies that for 
each solution x{t) of system (110.131) and to > 0 we have ||J/o) || < Bq for all t > to + T{a). 

The following condition is required: 

(Al) There exists a positive number Mg such that sup |l 5 (a::)|| < Mg. 

Theorem 10.1 Suppose that condition (Al) is fulfilled and there exists a Lyapunov function V{x) de¬ 
fined on R" such that V (x) has continuous first order partial derivatives. Additionally, assume that there 
exists a number B > 0 such that the following conditions are satisfied on the region ||x|| > B : 

(i) V (x) > a (||x||) , where a{r) is a continuous, increasing function defined for r > B which satisfies 
a{B) > 0 and a[r) —>■ oo as r ^ oo; 

(ii) (^) — ~^(ll^ll)’ where b[r) is an increasing function defined for r > B which satisfies 

b{B) > 0; 
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(iii) 


dV 

-g—{x) < c (||a:||), where c{r) is a function defined for r > B and there exists a positive number 

Mq such that 0 < c(r) < Mob{r) for all r > B. 


Then, for sufficiently small \p ,\, the solutions of system (110.1511 are uniformly ultimately bounded. 


Proof. Fix arbitrary numbers to > 0, a > 0 and a solution x{t) of system (110.131) . Take a number 
/3 satisfying 0 < /? < b{B). We consider system (110.1511 with a nonzero number ^ which satisfies the 
inequality 

MoMg {^~b(B)) ■ 

Our aim is to show the existence of numbers Bq > B and T{a) > 0, independent of Iq, such that if 
II 2 / 0 II < a, then ^ 0 , yo)|| < Bq for all t > to+T{a). 

Consider an arbitrary yo G R" such that Hyoll < ct- For the sake of brevity, let us denote y{t) = 
^Oj 2/o)- In the proof, both of the possibilities ||?/o|| < B and ||?/o|| > B will be considered. We start 
with the former. 

Let Mv = max V(x). Since a(r) —> 00 as r —^ 00, there exists a number Bn > B such that 
0'{Bo) > My ■ 

Now, suppose that there exists a moment si > tg such that || 2 /(si)|| > Bq. It is possible to hnd a 
moment S 2 satisfying tg < S 2 < si such that ||y(s 2 )|| = B and ||2/(t)|| > B for all t G [s 2 ,si]. 

Assumptions (ii) and (iii) imply for S 2 < t < si that 


dV(y(t)) dV, ^ 

-■ (/(2/(o) + l^9(x(t))) 

<-b(\\y(t)\\) + \t,\M,c(\\yit)\\) 

< (IpI MgMg - l)b(B) 

<-/ 3 , 


where “ • ” denotes the scalar product. 

The last inequality implies that F(y(si)) < V{y{s 2 )). On the other hand, by the help of assumption 
(i), we have V(y(s 2 )) < My < a{Bg) < F( 2 /(si)). This is a contradiction. Therefore, for all t > tg the 
inequality ||y(t)|| < Bg is valid. 

Next, we consider the possibility Hj/oH > B. Since the function V(x) is continuous and ||j/o|| < oc, one 
can find a number K(a) > 0 such that V(yg) < K(a). By means of condition (i) used together with the 
inequality Hj/qII > B, we have that K(a) > a(B). 

Assume that there exists a moment t > tg-\ - ^ ^ ^ ^^ such that ||2/(t)|| > B. 

If there exists ti G [tg,t] such that ||y(ti)|| < B, then by means of uniqueness of solutions, using a 
similar discussion to the case Hj/qII < B considered above, one can show that for all t >ti the inequality 
||?/(t)|| < Bg holds. On the other hand, if for all t G [tg,t\ the inequality ||y(t)|| > B is valid, then one 
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can verify that the inequality 


V{yit)) < V{yo)- P{t-to) 

holds. Under the circumstances we attain that 


a{B) < V{yo) - l3{t - to) < K{a) - f3{t - to) < a{B). 


This is a contradiction. Hence, for all t> to + T{a), where T{q) = ^ have ||y(t)|| < Bo- 


Consequently, the solutions of system (110.151) are uniformly ultimately bounded. □ 

Next, we shall verify the conditions of Theorem llO.ll for the Lorenz model. Let us consider the system 
(12.31) with the parameters ct > 0, 0 < r < ^/2 — 1, 6 > 0, and take into account the Lyapunov function 


V (u) — —Ui + U 2 + Uo, 
a 


where u = {ui,U 2 ,Uo) € 
1 


Set 7 i = min j^l, —j and define the function a(r) through the formula a{r) = In that case, the 
relation V{u) > 71 (rtf + + u\) = a (||u||) holds. On the other hand, one can verify that 


= —uiu'i + 2U2U2 + 2 uou'^ 

2 

= + M2)] + 2m2(—U1U3 + rui — U2) + 2 uo{uiU 2 — buo) 

a 

= 2 (r + 1 )miM2 — 2u\ — 2u\ — 2tbu\. 


Now, let 72 = min {l, 2 — (r + 1)^, 25} . Making use of the identity 

2 (r + 1)miM2 = u \ + (T + 1)^M2 - [ui - (r + 1)m2]^ 


we attain the inequality 

= -[ui-ir + 1 )m 2 ]^ - m^ - [2 - (r + 1 )^] uj - 2bul 
< —uf — [2 — (r + 1 )^] M 2 — 26 M 3 

where the function b(r) is defined through the formula b(r) = 72 r^. The last inequality validates the 
condition (ii) of Theorem 1 10. II 
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Furthermore, one can obtain that 


= ‘ 2 \l^ul+ul + ul < c(||u||), 

r 11 273 

where c(r) = 273 r and 73 = max < 1, — > . If we take Mq = -, then the inequality c(r) < M^bir) holds 

I cr J 72 

for all r > 1. Consequently, for i? = 1, the conditions of Theorem 110.11 are satisfied for system (12.311 with 
the coefficients ct > 0 , 0 < r < — 1 and 6 > 0 . 

In the next section, we will continue with the extension of sensitivity, which can be considered as the 
unique ingredient of chaos for a set of bounded solutions gsiEniEni- 

10.2 Unpredictability analysis 

Extension of the sensitivity feature through system (110.151) will be handled in the present part. We shall 
begin with the meaning of the aforementioned property for systems (110.131) and (110.151) . The main result 
will be stated in Theorem 110.21 

System (110.131) is called sensitive if there exist positive numbers eo and A such that for an arbitrary 
positive number Sq and for each chaotic solution x{t) of system (110.131) . there exist a chaotic solution x{t) 
of the same system and an interval J C [0, 00 ), with a length no less than A, such that ||a;(0) — x(0)|| < Jq 
and \\x{t) — x{t)\\ > cq for all t € J. 

Our main assumption is the existence of a bounded positively invariant set for system (110.151) . 
The existence of such an invariant set can be shown, for example, by using Theorem ll0.il 

We say that system (110.151) is sensitive if there exist positive numbers ei and A such that for an 
arbitrary positive number 61 , each i/q G and a chaotic solution x{t) of (110.131) . there exist yi G 
a chaotic solution x{t) of (|10.13|) and an interval C [0,oo), with a length no less than A, such that 
llyo - yill < <^1 and \\(t>x{t){t,0,yo) - 0 , 2 / 1 )II > ei for all t S J^. 

The following assumptions are needed: 

(A2) There exists a positive number Mp such that sup ||F(a:)|| < Mp; 

(A3) There exists a positive number Lf such that \\f{yi) — /( 2 / 2 )|| < Lf II 2/1 — 2 / 2 II for all 2 / 1 , 2/2 S M"; 

(A4) There exists a positive number Lg such that ||(/(a;i) — ( 7 (a; 2 )|| > Lg ||a;i — X 2 \\ for all a;i,a ;2 G K"*. 

In the next theorem, the extension of sensitivity from system (110.131) to system (110.151) is considered. 

Theorem 10.2 Suppose that conditions (Al) — (A4) hold. If system (110.131) is sensitive, then the same 
is true for system (110.151) . 

Proof. Fix an arbitrary positive number i5i, 2/0 S and a chaotic solution x{t) of (110.131) . Since system 
(I10.13|) is sensitive, one can find eo > 0 and A > 0 such that for arbitrary > 0 both of the inequalities 
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||a;(0) — s(0)|| < ^0 and ||a;(i) — s(<)|| > cq, t G J, hold for some chaotic solution x{t) of (110.131) and for 
some interval J C [0, oo), whose length is not less than A. 

Take an arbitrary yi G such that ||?/o ~ 2/i|| < For the sake of brevity, let us denote ?/(f) = 
</>a:(t)(i,0,?/o) and y{t) = (j>^(^t){t,0,yi). 

It is worth noting that there exist positive numbers Kq and Hq such that ||?/(t)|| , ||y(t)|| < Ko for all 

t > 0 and sup ||a;(t)|| < Hq for each chaotic solution x{t) of system (110.131) . 
t>o 

Our aim is to determine positive numbers ei, A and an interval C [0, oo) with length A such that 
the inequality \\y{t) — y{t)\\ > ei holds for all t G J^. 

Since the derivative of each chaotic solution x{t) of (I10.13|) lies inside the tube with radius Mp, the 
collection of chaotic solutions of system (110.131) is an equicontinuous family on [0, oo). Suppose that 
g{x) = {gi{x)^ g 2 {x),..., gn{x)), where each 1 < j < n, is a real valued function. Making use of 
the uniform continuity of the function g : R™ x R™ —?► R", defined as g(^ 1 , 1 / 2 ) = g{vi) — g{v 2 ), on 
the compact region = {(I'l, 1 ^ 2 ) G x : ||i^i|| < Hq, \\v 2 \\ A .^o} together with the equicontinuity 
of the collection of chaotic solutions of (I1Q.13L one can verify that the collection ^ consisting of the 
functions of the form gj{xi{t)) — gj{x 2 {t)),l < j < n, where a:i(t) and X 2 {t) are chaotic solutions of 
system (ll().13L is an equicontinuous family on [0, 00 ). 

According to the equicontinuity of the family ^, one can find a positive number r < A, which is 
independent of x{t) and x{t), such that for any ti,t 2 G [0, 00 ) with \ti — t 2 \ < r, the inequality 

\{9j (xih)) - gj (xih))) - {gj {x{t 2 )) - gj (x(t 2 )))| < (10.16) 

holds for all 1 < j < n. 

Condition (A4) implies that ||g(x(t)) — g{x(t))\\ > Lg \\x{t) — x(t )||, t G J. Therefore, for each t G J, 
there exists an integer jqi 1 < jo < which possibly depends on t, such that 

\9joix{t)) - 9joix{t))\ > ^ \\x{t) - S(t)|| . 

Otherwise, if there exists s G J such that for all 1 < j < n the inequality 

li'j (a;(s))-5j{^{s))l < ^ lk(s)-^(s)|| 


holds, then one encounters with a contradiction since 

n 

||g(x(s)) - g(s{s))|| < l5i(a;(s)) - gj(x(s))| < Lg ||a:(s) - x(s)||. 

i=i 

Denote by Sq the midpoint of the interval J, and let 9 = sq — r/2. There exists an integer Jq, 
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1 ^ jo < n, such that 

l3jo(a^(so)) -5io(^(so))| > ^ ||a;(so) -S(so)|| > (10.17) 

On the other hand, making use of the inequality (110.1611 it can be verified for alH G [0,0 + r] that 


\93o {x{so)) - gj„ (a;(so))| - \9jo {x{t)) - gj^ (a;(t))| 

< \ i9jo (a^(0) - 9jo - {9io (a^(so)) - 9jo (^(so)))| 

Lg€0 
2 n ■ 

Therefore, by means of (I10.17L we have that the inequality 

\9jo - gj^ (x(t))| > \gjg {x{so)) - gjo (x(so))| - (10.18) 

is valid for t G [0,0 + r]. 

One can find numbers si, S 2 ) ■ • ■ j G [0,0 + r] such that 

_ _ 

/ b( 2 ;(s)) - 0 (a;(s))] ds = (^T [gi(a::(si)) - gi(x(si))] ,t [g 2 (a;(s 2 )) - g 2 (a;(s 2 ))] , 

^ & 

• ■ • , T [5n(x(s„)) - 5n(^(Sn))] ) ■ 


By using the inequality (110.181) . we attain that 


10 


[ 5 (a;(s)) - 5(a;(s))] ds 


> T \gjo{x{sj^)) - 5jo(®(sio))l > 


xLgCo 

2 n 


The relation 


y(0 - y{t) = {y{d) - y{0)) + [ [fivis)) - f{y{s))] ds+ f g.[gix{s)) - g{x{s))]ds, t G [0,0 + r] 

Js Jo 


yields 


|ly(0 + T) -y(0 + T)|| > 


rd+r 


[0(a;(s)) -5(2;(s))]ds 


- \\yi 0 ) - yi 0 )\\ - f Lf b(s) - y(s)|| ds 

Js 

- r\f\\y{s)-y{s)\\ds. 

Js 


2 n 
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The last inequality implies that 


max \\y{t) - y{t)\\ > \\y{e + t) - y{e + t)|| 

te[e,e+r] 


ImI ThgCp 
2n 


(l + rLf) m&x \\y{t)-y{t)\\. 

t^\u ^u-\-t\ 


Therefore, max \\y(t) — y{t)\\ > 
tgfe.e+r] 


ImI TLgep 
2n(2 + rLf) 


Suppose that max \\y{t) — y(t)|| = \\y{0 ~ ?/(C)ll for some ^ S [0,9 + t]. Define 

tGl0,e+T] 


A = min 


(t _ \fi\ rTg£o _ 

I2’ 8niKoLf + Mg\y\){2 + TLf) 


and let 


= 


ii^<e + T/2 

^ - A, if e > 0 + r/2 


For t G + A], by favour of the equation 


yit) - y{t) = (yiO - yiO) + J Ifiyis)) - f{yis))]ds + J y[g{x{s)) - g{x{s))]ds, 


one can obtain that 


lll/W-yWII> 112/(0-^011- 

- \y\ Il5(a^(s)) -5 (^(s))Ns 

\y\TLgeo 


Lf\\y{s)-y{s)\\ds 


> 


> 


2n(2 + rLf) 
\y\rLgeo 
4n(2 + rLf) 


-2A {KoLf + Mg\g\) 


The length of the interval = [0^, 0^ + A] does not depend on x{t), x{t), and for t G the inequality 

\y\ TLgCO 


\\y(t) — y{t)\\ > ei holds, where ei = 


4n(2 + rLf) 


. Consequently, system (110.1511 is sensitive. □ 


10.3 Existence of unstable periodic motions 

Assume that system (110.131) admits a period-doubling cascade. That is, there exists an equation 


x' = G{x,\), (10.19) 

where A is a parameter and the function G : K.'" x K —?► R™ is such that for some finite number Aoo, 
G{x, Aoo) is equal to the function F{x) in the right hand side of system (ll().13|l . 
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System (110.131) is said to admit a period-doubling cascade pTl 15^ [S5l ITT] if there exists a sequence of 
period-doubling bifurcation values satisfying Xj —^ Aoo as j —^ oo such that as the parameter A 

increases or decreases through Xj system (110.191) undergoes a period-doubling bifurcation for each j G N. 
As a consequence, at the parameter value A = Aoo, there exist infinitely many unstable periodic solutions 
of system (110.191) , and hence of system (110.131) , all lying in a bounded region. 

Now, let us introduce the following definition |7D]. We say that the solutions of the non-autonomous 
system (110.151) . with a fixed x(t), are ultimately bounded if there exists a number B > 0 such that 
for every solution y{t), y{to) = yo, of system (110.151) . there exists a positive number R such that the 
inequality ||y(t)|| < B holds for all t > to + R. 

We say that system (110.151) replicates the period-doubling cascade of system (110.131) if for each periodic 
solution x{t) of (110.131) . system (I10.15|) admits a periodic solution with the same period. 

The following condition is required in the next theorem, which can be verified by using Theorem 15.8 

m- 

(A5) Solutions of system (I10.15|) are ultimately bounded by a bound common for all x{t). 

Theorem 10.3 If conditions (Al) — (A5) hold, then system il().15\} replicates the period-doubling cascade 
of system il(J.lS\) . 

It is worth noting that the instability of the inhnite number of periodic solutions of system (I10.15P is 
ensured by Theorem 110.21 
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